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ABSTRACT 


In this paper, a high accuracy linear positioning system with a linear force actuator and magnetic 
levitation is proposed. By locating a permanently magnetized rod inside a current-carrying solenoid , the 
axial force is achieved by the boundary effect of magnet poles and utilized to power the linear motion, while 
the force for levitation is governed by Ampere's Law supplied with the same solenoid. With the levitation 
in a radial direction, there is hardly any friction between the rod and the solenoid. The high speed motion can 
hence be achieved. Besides, the axial force acting on the rod is a smooth function of rod position, so the 
system can provide nanometer resolution linear positioning to the molecule size. Since the force-position 
relation is highly nonlinear, and the mathematical model is derived according to some assumptions, such as 
the equivalent solenoid of the permanently magnetized rod, so there exists unknown dynamics in practical 
application. Thus "robustness" is an important issue in controller design. Meanwhile the load effect reacts 
directly on the servo system without transmission elements, so the capability of "disturbance rejection" is 
also required. With the above consideration, a time-delay control scheme is chosen and applied. By comparing 
the input-output relation and the mathematical model, the time-delay controller calculates an estimation of 
unmodeled dynamics and disturbances and then composes the desired compensation into the system. 
Effectiveness of the linear positioning system and control scheme are illustrated with simulation results. 


INTRODUCTION 


Interest in research on large-gap magnetic suspension systems began in the early 1960's. The 
principal goal was the elimination of aerodynamic support interference in wind tunnel testing. To early 
1970 s the interest extended to small-gap ones. The first system developed was the Annular Momentum 
Control Device (AMCD) with applications to the stabilization and control of spacecraft [1], This research 
was continued with the Annular Suspension and Pointing System (ASPS) which provides orientation, 
mechanical isolation, and fine pointing of space experiments[2][3]. For decades, Magnetic suspension 
technologies (MST) have demonstrated their capabilities in many fields, from industrial compressors, 
high-speed milling and grinding spindles, magnetically levitated trains, control wheel suspension for 
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spacecraft to rocket propulsion turbomachinery. Important features of the magnetic suspensions and 
actuator systems are. 

(1) Versatility of the Electromagnetic Forces 

The physical force of a magnetic circuit to a high-permeable armature is called the Maxwell-force. 
Contrary to this commonly used force, the reaction force of a conductor carrying a current in a 
magnetic field is called the Lorentz-force. By successful integration of these physical effects, the 
constructed electromagnetic subsystem can be utilized as a rotary motor, linear actuator, radial 
bearing, thrust bearing, etc. 

(2) Molecule-size Resolution 

One problem of electric motors is the ripple of motion at low-speed operating regions due to the 
finite pole effect. The rotor always rests at the finite circumference positions which have the 
minimum magnetic flux (potential energy). Thus there are inherited limitations for resolution of 
control. The non-pole magnetic field provided by a coil, on the other hand, sets no resolution 
limitation. The resolution limit, in turn, is set by sensors, instrumentation and control strategies. 
Magnetic suspension systems provide a promising approach for achieving positioning with 

nanometer resolution. ... 

In this paper, a linear positioning system with a linear force actuator and magnetic levitation is to be 

designed. By locating a permanently magnetized rod inside a current-carrying solenoid, the axial force is 
achieved by the boundary effect of magnet poles and utilized to power the linear motion, while the force for 
levitation ispiovidedby the magnetic bearing and governedby the maximum linkageprinciple.With the levitation 
in a radial direction, there is no friction between the rod and solenoid. The demand of high-speed motion can 
hence be achieved . Under the proposed arrangement, the axial force acting on the rod is a smooth function of 
rod position, so the system can provide nanometer-resolution linear positioning to the molecule size. It is 
known that an effective control system is the key condition for successful implementation of 
high-performance magnetic levitated positioning systems. Major issues for design of such control systems 

are: 


(1) Nonlinearity . . . 

By assuming that the complete energy of the magnetic field is concentrated within the air gap. 1 he 

basic mathematical models of active magnetic bearings are obtained from Maxwell's laws. The 

input-output relations are highly nonlinear despite the variables defined. 

(2) Unmodelled Dynamics 

Secondary effects such as copper resistance, stray fields and saturation are neglected. 

(3) Disturbance Rejection ... 

Because the load effect reacts directly on the servo system without the transmission elements, the 

capability of "disturbance rejection" is also required. 

With the above considerations, a time-delay control scheme with the properties of "robustness" and 
"disturbance rejection" is utilized[4][5]. By comparing the input-output relation and the mathematical 
model, the time-delay controller calculates an estimation of unmodeled dynamics and disturbances and then 
composes the desired compensation into the system. Effectiveness of the linear positioning and control 
systems is illustrated by numerical simulation results. 
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SYSTEM DYNAMICS 


System Configuration 


The configuration of the proposed magnetic levitated linear positioning system is shown in Figure 1. To 
achieve the function of levitation, the current in the solenoid must be 'kept in the direction that can maintain 
the stability of radial motion. Under such condition, the magnetic force in the axial direction tends to push the 
rod away from the center of the solenoid. Hence a spring is required to supply the force in the opposite 
direction. Also, the spring must be precompressed to avoid an uncontrollable equilibrium point. The additional 
magnetic bearing system is used to keep the moving part balanced in the axial direction. With a biased current 
fed to the solenoid, the magnetic force ( F=Ixfj )in the radial direction is utilized to suspend the moving 
part, while, with the controlled current, the axial motion is governed by the force caused by the non-uniform 
magnetic field in the boundary. 


Permanently 



Figure 1 The configuration of magnetic levitated linear positioning system 


Dynamics in Axial Motion 


To derive the dynamics, a solenoid is employed to produce the equivalent magnetic field B of the 
permanently magnetized rod and thus simplify the calculation of the force-position relationship. With such an 
arrangement, the solenoid-rod configuration can be approximated as depicted in Figure 2. The force acting 
on the inner coil due to the current in the outer coil can then be found with Ampere's Law. 
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Figure 2 Equivalent configuration 


Consider two current-carrying elements as shown in Figure 3. The force between the two 
current-carrying elements is[6] 


d¥ - 


iao/ 2/1 

47tr 2 


ai 


x (a 2 x r)dl 2 dl\ 


where dF = force on element 1 due to current in element 2, N 
|io = permeability of air, Hm 1 

dl\,dl 2 = lengths of current-carriyng elements 1 and 2, respectively, m 
I \, 1 2 = currents in elements 1 and 2, respectively, A 
r = distance between elements, m 
a, = unit vector in direction of current in element 1, dimensionless 
a 2 = unit vector in direction of current in element 2, dimensionless 
r = unit vector in radial direction( from element 2 to 1), dimensionless 
Applying the force-position relation to the equivalent model, it can be found that the force acting on the 
inner coil with radius r 2 at position x 2 due to the current-carrying outer coil with radius r , at position x, is 
given by the following equation: 


dF = Hst! — . [cose J - Sin02*] 

4tc[(xi -X2) 2 +(risin0i - /• 2 sin0 2 ) 2 + (/"iCOsGi -r 2 cos0 2 ) 2 J 

x{[cos0j -sin0i*] x[(x 2 -xi)7 +(r 2 sin0 2 -nsinOi)/ +(r 2 cos0 2 -ricosGi)^ } 


Hence the force in the axial direction can be expressed as 

-po/i/ 2 rir 2 (x 2 -Xi)(cos0iCQS0 2 + sinOisin0 2 ) dQ\dQ 2 

4ti[(xi -x 2 ) 2 + (risin0i -r 2 sin0 2 ) 2 +(ricos0i -r 2 cos0 2 ) 2 ] 

with the numerical calculation, the force-position relation of the rod, with 5000 Gauss of flux density at 
one end on the axis shown in Figure 4, where r, is 1.1 cm, r 2 is 1.0 cm, the length of the rod is 1.0 cm, 
and the length of the solenoid is 10 cm. The result indicates that the solenoid-rod configuration 
operates as a current-controlled stiffness nonlinear spring, so the model of axial motion can be simplified 
and shown in Figure 5, where K, is the preload spring, and is the equivalent model of the solenoid-rod 
configuration. 
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Figure 3 Geometry of short current-carrying elements for finding force between them. 
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Figure 4 Force-position relation in axial motion. 
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The dynamic equation in axial motion can then be modeled as 

Xi = x 2 

x 2 =-jj[Ki(xi +x p )+K 2 (xi)(i + ib )] 


where 


x„x 2 

K 2 (x) 


/ 

m 


position and velocity of the rod respectively, cm, cm/sec 
stiffness of the spring, N/cm 

current controlled stiffness of solenoid-rod configuration, N/A 

pre-compressed length of the spring, cm 

biased current for levitation, A 

controlled current, A 

mass of the magnetic rod, kg 


Dynamics of Radial Motion 


To verify the function of levitation of the proposed mechanism, the dynamics of radial motion is taken 
into consideration. Imagine a small deviation of the rod position in radial direction^, then the radial force 
can be modeled as 

Un/i/2nr2k2sin02(cos0icos02 +sin0isin02) + (y2sin0] -ri)sin02] „ .q 

ifp — — — “ ' 13/2 1 y 

y 47i[(xj - x 2 ) 2 +(>2 +r,sin0i -r 2 sin0 2 ) 2 +(r,cos0i -r 2 cos02) 2 J 

From the calculation, the result is shown in Figure 6 . It is clear that, with positive current, the force 
tends to force the rod staying in the central axial position. Any deviation in radial direction will cause a 
force in the opposite direction to push it back. Hence the function of levitating is achieved. In this paper, 
the motion in the radial direction is not controlled. 

Since the system is passive and motion in the radial direction is confined, the open loop in the radial 
direction poses no stability problem. 
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CONTROLLER DESIGN 


Time Delay Controller 


Since the magnetic levitated positioning system is subjected to the load effect directly, the controller 
must be good at disturbance rejection. Besides, the force-position relation is highly nonlinear. Hence the 
Time-Delay Controller (TDC) design method proposed by Youcef-Toumi [3] is applied to control the rod 
position. 

Consider the axial motion to be controlled as 

x = f(x, /) + h(x, t) + B(x, /)u + d(/) 


where 


x = an n x 1 plant state vector, 

f (*,0 — an n x 1 nonlinear vector represents the known dynamics, 
h(x,/) = an n x 1 nonlinear vector represents the unknown dynamics, 
®(x,0 = an n x r control distribution matrix with rank r, 
d(7) = an // x 1 unknown disturbance function. 


Define the linear time-invariant reference model as 
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and the error between the model and the plant as 


e = x„, - x. 

The objective of the control scheme is to make the error behave as the desired dynamics 

e = A e e = ( A m + K )c 

Combined with the plant dynamics and reference model dynamics, the error dynamics can be 
formulated as 

e = A e e + [-f - h - d +A m x +B„,r -Bu -Ke] 

To satisfy the desired error dynamics, the control input u must be functioned to make 

[-f - h - d + A„,x + B m r - Bu - Ke] = 0 

Observe the above equation; only h and d are unknown. They can be estimated from the measured 
information at time t —L, with 

h(x, 0 + d(0 = h(x, t-L) + d(t-L) = x(7 -L)-f(x,t-L)- B(x, t - L)u(t - L) 

Then the TDC control law is 

u(/> = B + (/)[f (x, t-L)-x(t-L) + B(x, t - L)u(t -L)- f (/) + A m x(/> + B„,r(0 - Ke(0] 

Rewrite the dynamics equation of the magnetic levitated linear positioning system as 

Xi= x 2 

x 2 = -±j{[Ki(x l + x p )] + [/: 2 (xi>- (/ + i b ) + mbi \ } + bi 


and let 


fX*) = 


*2 

h(x) = 

0 

A 

and B = 

’ 0 
A 

—K\ • (xi +Xp)//w 


-K 2 (x i) • (/ + i b )/m - bi 


b 


then by choosing the reference model as 


Xi 


0 

1 

Xi 

+ 

0 

_ X 2 


-100 

- 

-20 

X2 _ 


100 


and K = 0, then the control force is 
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/(/) = b 1 |aTj [*i (t - L) - x\ (t)\/m - X 2 O - L) + b(t -L)i(t ~L)~ 100xi(/)-20x 2 (/)+ 1 OOr(t) J 

The simulation results are shown in Figure 7 and Figure 8. In the simulation, the control force is limited 
with a saturation bound [-3.5, 10] and b is chosen to be a constant 10. Under the case in Figure 7, the 
linear positioning system is subjected to a step disturbance in the axial direction at time t = 0.5 sec, whose 
amplitude is -0. 1 . For the case in Figure 8, the system is utilized to move a load with mass equals to 5 m, 
the mass of the moving part. 






Figure 7 Simulation result with Time Delay Controller. 
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Figure 8 Simulation result of TDC 


Time Delay Controller with Mean Value Theorem 


The method described in the previous section requires the measurement of xi- It means that the designer 
must take care of either the noise induced by numerical derivativeness or the setting up of the 
accelerometer. An alternative scheme is introduced by Youcef-Toumi[4] with the application of 
the mean- value theorem as follows: 

Define the unmodeled dynamics - h + d; the error dynamics can be written as 

e = A e e + [-f -'F +A m x +B„,r -Bu -Ke] 

Integrate the equation at both sides to get 

e(r) = 0(7, t 0 )e (to) + f <*>(', t) [A m x(x) + B,„r(x) - fXx) - T(x) - Bu(x) - Ke(x)]^x Vr > to 

^ h 
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where 


0(7, x) = e A,(r t) 

Then 

j' o 0(7, x) v F(x)Jx = 0(7, t 0 )t(to) - e(0 + 0(7, x) [A m x(x) + B m r(x) - f(x) - Bu(x) - Ke(x)]rfx 

Utilizing the mean-value theorem approximation, the above equation can be approximated by 
[j' o 0(/, x) dx 'F(t) = 0(7, / 0 )e(/ 0 ) - e(0 + 0(7, x) [A m x(x) + B m r(x) - f(x) - Bu(x) - Ke(x)]afc 

Thus the control force can be found as 

u = B + [A m x +B m r -f -if -Ke] 

With the same choice of reference model in the previous section, the simulation result with step 
disturbance in axial force is given in Figure 9. Figure 10 demonstrates the result of load effect. 



time (sec) time (j ac ) 

Figure 9 Simulation result with modified Time Delay Controller (I) 
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Figure 10 Simulation result with modified Time Delay Controller (II) 


Discussions 


From the simulation results, the property of "disturbance rejection" of Time Delay Controllers is 
demonstrated. Due to the saturation constraint in the control force, the controller can perform well with 
limited external load. With the mean-value approximation, the performance is improved. But it requires 
a higher sampling rate to gather more information in the delay time. 

During the design stage, the choice of parameter b plays an important role. In general, if we know 
more about the nonlinear dynamics of the system to be controlled, we can have a better chojce of b. Under 
the situation that no information is given for the unknown dynamics, we usually choose b as a constant 
value. From simulation results shown in Figure 1 1, we conclude that the larger b provides more guaranteeon 

the "disturbance rejection", but has poorer performance. But, the larger b drives the control force to the 

saturation bound easily and consequently the system fails to work. 
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Figure 1 1 Response under different choices of b . 


CONCLUSION 


Applying the equivalent magnetic flux density between the solenoid and permanently magnetized rod, a 
method of modeling the force-position relation in the solenoid-rod configuration is proposed. Based on the 
model, a "current controlled stiffness" spring serves as a low friction linear positioning system. Two kinds 
of Time Delay Control schemes are applied to control the system. The Time Delay Controller is a fast 
adaptive control scheme without identification of the dynamics of the plant to be controlled, so it can be 
easily implemented. From the simulation results, the property of disturbance rejection is demonstrated. 
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